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Abstract 

A technique  is  developed  whereby  one  can  syn- 
thesize a causal,'  linear  tit e- invariant  estimator 
that  is  optimal  for  restricted  types  of  nonstation- 
ary processes.  The  technique  is  applicable  to 
linear,  t irce- invariant  systems  (driven  by  nonsta- 
tionary state  noisp)  for  which  scalar  observations 
ere  cade  in  the  presence  of  additive  nonst at ionax y 
noise.  Two-dimensional  Fourier  transforms  are  used 
to  obtain  an  expression  for  the  estimator's  mean 
square  error.  It  is  assumed  that  it  is  desirable 
to  minimize  the  time  integtal  of  this  expression. 
The  calculation  of  this  integral  results  in  an 
expression  which  car.  be  minimized  by  selecting  an 
estimator  depending  in  a prescribed  way  on  the 
two-dimensional  Fourier  transforms  of  the  state  and 
observation  noise.  The  resulting  es tier. tor  is 
causal,  linear,  and  time  invariant.  It  is  similar 
in  some  respects  to  the  Wiener  filter  that  can  be 
derived  under  the  assum.pt icr.s  of  stationary  state 
and  observation  -noise  processes.  The  esiir.acor'* 
usefulness  is  limited  by  the  requirement  that  the 
observations  be  scalar,  and  the  r.onstatior.ary  pro- 
cesses have  Fourier  transformable  autocorrelation 
functions. 1 


I.  Introduction 

The  determination  of  optimal  linear  time 
invariant  estimators  that  can  be  used  to  process 
the  output  of  a linear  tine  invariant  system  driven 
by  stationary  noise  and  observed  in  the  presence  of 
additive  stationary  noise,  was  the  subject  of 
Wiener's  work  [l].  This  early  work  was  followed  by- 
Kalman's  [2]  determination  of  a method  for  finding 
optimal  estimates  given  certain  types  of  nonsta- 
tionary disturbances  and  corruptions.  /.  major  dis- 
tinction between  these  two  results  was  the  methods 
by  which  they  were  obtained.  Wiener's  method  was 
frequency  domain  based,  while  Kalman's  approach 
relied  on  time  domain  techniques.  A survey  of  per- 
tinent literature  indicates  chat  until  Kalman's 
results  were  published,  the  major  work  on  nonsta- 
tionary processes  was  confined  to  looking  at  slight 
modifications  of  the  stationary  problem. 2 When 
Kalman's  results  were  announced,  there  was  a sodden 
transition  from  frequency  domain  methods  and  little 
attention  has  been  paid  to  frequency  domain  techni- 
ques in  estimation  since  that  rLrnc. 

In  many  applications,  the  Kalman  estimator 
results  in  exhorblta.it  computational  requirements. 


This  is  due  mostly  to  the  fact  that  its  time  vary- 
ing covariance*  matrix  must  be  calculated  before 
the  estimator  gains  can  be  computed.  In  an  effort 
to  cake  the  estimator  less  complex,  stost  designers 
use  the  steady  state  gain  values  or  schedule  sets 
of  gain  values  that  are  to  be  used  at  different 
points  it;  the  processing  algorithm.  This  results 
in  a simpler  estimator  that  performs  suoopilmal ly 
during  periods  when  the  process  being  observed  is 
nonstationary . The  estimator  will  be  optimal  when 
the  process  becomes  stationary  provided  the  steady 
state  value  of  the  estimation,  error  covariance  is 
used  in  the  gain  computation. 

Use  of  the  stationary  gain  values  constitutes 
selection  of  the  Wiener  estimator  as  a special 
case  cf  the  Kair.an  estimator.  It  is  this  proce- 
dure fox  selecting  a subopt ir.al  estimator  that 
prompts  one  to  ask  if  there  might  oe  a causal, 
linear,  time-invar iar.r  estimator  that  is  optimal 
in  some  restricted  sense  for  seme  class  of  non- 
stat'or.ary  processes. 

To  arrive  at  such  an  estimator,  it  seems  rea- 
sonable to  consider  the  case  most  likely  to  yield 
such  a result.  Namely,  the  system  or  state  model 
should  be  linear  and  tine  invariant.  It  should  be 
driven  by  zero  mean  r.onstatior.ary  r.c ise  and 
observed  in  the  presence  of  additive  nonstat icr.ary 
noise.  Optimality,  for  the  sake  of  obtaining 
‘results  comparable  to  Wiener  and  Kalman,  should  be 
defined  in  the  sense  that  the  mean-square  estima- 
tion error  or  some  function  of  it  is  minimized. 

If  cne  returns  to  frequency  domain  considerations, 
it  might  be  possible  tc  arrive  at  an  estimator  via 
techniques  paralleling  those  of  Wiener,  provided 
Fourier  transforms  of  the  noise  proc  .sses  can  be 
defined.  Tr.is  leads  to  the  use  of  two-dimensional 
Fourier  transforms  to  obtain  the  average  energy 
spectral  densities  of  the  state  and  reservation 
ncisc-  processes.  An  expression  is  arrived  at  fer 
the  rear.-square  estimation  error  and  it  is  a* fumed 
that  the  optimal  estimate  should  minimize  the 
doubly  infinite  time  integral  of  this  error.  The 
expression  that  results  from  doing  this  is  used  to 
determine  th  structure  of  the  optimal  estimator. 
The  resultant  estimator  Is  similar  in  sore 
respects  to  the  Wiener  estimator.  The  estimator 's 
usefulness  is  limited  by  the  requirement  that  the 
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observat  lours  be  scalar  and  the  noise  processes 
have  Fourier  transformable  autocorrelation  func- 
tions thnl  can  successfully  undergo  spectral  fac- 
torization. 

11.  System  Assumptions  and  Outline  of  Approach 

The  solution  of  the  problem  posed  in  the 
introduction  is  developed  in  four  steps.  These 
are: 

1.  Consider  a linear,  t ime- Invariant  system 
driven  by  nonstationary  noise  and  observed  in 
the  presence  of  additive  nonstationnry  noise. 
Assume  the  observations  are  scalar. 

2.  Formulate  an  expression  for  the  mean  square 
estimation  error  of  a linear  time-invariant 
estimator  with  impulse  response  h(T),  that 
might  be  used  to  obtain  an  estimate  of  the 
noise-free  observation. 

3.  The  expression  for  the  mean  square  estimation 
error  obtained  in  step  2 will  in  general  be 
time  dependent  due  to  the  noise  being  nonsta- 
tionary. Provided  that  one  is  willing  to 
accept  a filter  that  minimizes  this  expression 
in  an  average  sense,  i.e.,  over  a doubly 
infinite  time  interval,  one  can  arrive  at  a 
simpler  expression  for  the  integral  of  the 
mean  square  error  that  does  not  depend  on 

t iae . 

4.  Finally,  if  one  restricts  the  type  of  nonsta- 
tionary noise  processes  sufficiently,  the 
expression  obtained  in  step  3 can  be  minimized 
by  a properly  designed  causal,  linear,  time- 
invariant  estimator  that  minimizes  the  inte- 
gral of  the  .mean  snuarp  estimation  error. 


It  .is  desired  to  determine,  for  every  t,  an 
estimate  D(t),  for  w( t)  which  minimizes 

F.[e)(t)  - £>(t)]?  (3-3) 

or  some  function  of  this  expression  by  processing 
only  the  "present  and  past"  values  of  the  observed 
data: 

y(T)  - s(t)  + n(T)  - *>  £ T £ t 
It  is  further  required  that 

CD 

Hj(t)  - / h(T)y(t-T)  d-.:  h(T)  - O.t  < 0 (3.4) 


In  preparation  for  arriving  at  an  expression 
for  the  estls.aior's  mean  square  estimation  error 
let  us  introduce  the  tvo-diciens ional  Fourier 
tr  ansfornu 

Definition:  The  two-dimensional  Fourier  transform 

of  a function  fCt^t^)  is  given  by:2 

T?  -j^t-het) 

FCJUpjup  * j \ f(trt2)e  dtl  dt2 


its  inverse  is  defined  by 


1 f C * 1 1 

f(t  ,t  ) l—r  / / F(ju ..Jo ,)e 


j(u.t,  + ce,t-) 


d. do.' j 


The  following  sections  will  expand  and  elabo- 
rate on  the  four  steps.  A final  section  will  dis- 
cuss the  implications  of  the  final  result. 


As  in  the  case  of  one-dimensional  transforms, 
these  two  expressions  constitute  a Fourier  trans- 
form pair  that  we  will  represent  symbolically  as 


III.  System  and  Noise  Models 

Proceeding  with  the  problem  solution,  consider 
the  error  propagation.  The  system  being  considered 
is  linear  and  time  invariant.  As  such,  it  is 
described  by  the  following  linear,  time-invariant 
vector  matrix  differential  equations: 

x(t)  - A_x(t)  + Bii(t)  (3.1) 

y (t)  - Cx(t)  + Dv(t)  - s(t)  + n(t) 


_x(t):  is  a vector  representing  the  stare  of 

the  system 

y(t):  Is  a vector  representing  an  observation 

of  the  state  of  the  system  (scalar) 


f (t.  ,t7)  I — ► F(Ju>.  ,Je!  ) 


It  should  be  noted  that  if  f(tj,t2)  is  such 


® OO 

J j | f(tx.t2)  | dtxdt2 


then  FfjUj.ju.’j)  exists  for  every  srd  11  ls 

bounced . 


u(t),v(t):  are  random  processes 


The  problem  Is  to  estimate  any  time-invariant  (but 
not  necessarily  causal)  linear  functional  w(t)  of 
the  signal  s(t)  defined  by  the  following  infinite 
convolution  integral: 


w(t)  - J 1(A)  s ( t - A)  dA 


JFor  a description  of  two-dimensional  Fourier 
transforms  and  their  application  to  nonstationary 
processes,  refer  to  Papoulis,  A.,  Probab 1 i i t y , 
Random  Variables,  and  Stochastic  Prorer ses. 
McGraw-Hill,  1565,  Chapter  12. 


IV 


The  Mean  Square  Error  Expression 


with 


Uc  can  compute  an  expression  for  r(t)  “ iu(c) 
£)(t).  To  do  this  recognize  that: 


u(t)  - J y(t-i)  h(T)  Jx 

m-eo 

and  for  the  present  let 

<a(t)  - J s(t-r)  1(t) 


dx 


(4.1) 


(4.2) 


then  by  subtraction 

e(t)  - <o(t)  - £(t) 
and  the  autocorrelation  of  e(t)  is 

Ree(tl.t2)  ■ EfeCt^eCt^  } - ^{[uCtjl-fiCtj)] 


[u(t,)-u(t2)] ) 


(4.3) 


tf(Wj)  * 

t 

W(u>2)  & 

being  the  standard  one-dimenslonal  Fourier  trans- 
forms of  U(tj)  ar.d  W(t2)  respectively. 

Taking  the  two-dimensional  Fourier  transform 
of  Ree(cj.t2)  we  have: 

ree<“i»“2)  " rSs(ll,i’tl!2) 

[i(j«2)-H(3«2)]  ♦rM(«l.«2) 
[hOUj)]  [h(Jl:2)J  (4.8) 


-JUjt, 


J V(t j)e  . 'dtj  - W(Jajj) 


-Js1 


(4.7) 


J W(t2)e  2 2dt2  - W(Ju2) 


Substitution,  expansion,  and  the  interchange  of 
expectation  and  integration  yields: 


If  we  now  take  the  inverse  transform  of  r 
we  have  an  expression  for  the  autocorrelation  of 
the  estimation  error: 


Ree<tl-t2)  * Rss<tl‘t2)*i(tl)*i(t2) 

- Rss(t1,t2)*i(t1)*h(t2) 

- ».,«*1-t2)*h<tl )*Kt2) 


od  at>  / 

/ / rss(“r“25  ' 

— ot  — ao  v 

5'KJjupJ  ^I(ju)2)-H(jm2)j 


+ S>s(t1,tj5*h(i1)*h(t2) 


+ Rnn(t1,t2)*h(t1)*h(t2)  (4.4) 


d terms  in  R an<^ 

and  denotes  convolutic 


Where  we  have  assumed 
R s(ci»c2^  are  zero  and  "**'  denotes  convolution  in 
tRe  time  domain.3  Assume  each  of  the  above  corre- 
lation functions  has  a two-dimensional  transform 
defined  as  follows: 


hI 


Rss(tl,t2)  ~ rss(“l’“2) 
*2 


“■1 

Rnn(tl,t2)  ~ rnn(ul,u2) 
*2 


(4.5) 


Kext  we  make  use  of  a property  from  £3}  page  442. 


R(tj,t2)*W(t1)*U(t2)  tZ  rfUj.cPjJWfUjJW^) 

*5 

2 (4.6) 


^“iV^V  r i 

* +rnn<“l-“2HR<JUPJ 


r t ^ ^i1  i+tL2t2^  I 

[H(Ju2)Je  jdujdu,  (4.9) 

Recall  that  we  seek  an  expression  for  the  mean 
square  estimation  error  for  each  point  in  time. 
Clearly,  if  t.  * tj  * t in  the  previous  expression 

we  have: 


Kee(t,t)  - E{e(t)e(t) ) 


(2x> 


C 09  1 

r J(u-Hj.)t 

[l(jL:2)-H(ju2)je 

r ir  i ) (|*)i+u)2)t 

[H(jUj)JILH(jv2)Je 


dUjduij 


(4.10) 


*1116  cross  terms  can  be  carried  through  the 
development  but  they  provide  no  additional  prob- 

lems and  therefore  have  teen  deleted. 


**v 


Note  that  this  expression  for  the  n 'in  square  esti- 
mation erroi  is  completely  penernl  except  for  the 
restriction  that  the  indicated  fir.nal  and  noise 
two-dimensional  transforms  exist.'  • 


If  we  ponder  this  expression  for  a moment,  we 
see  that  to  select  an  estimator  H(ju)  that  would 
minimize  it  for  each  value  of  t,  is  not  a simple 
task.  In  fact,  it  is  not  obvious  that  such  an 
estimator  even  exists. 


It  i s Interesting  to  note  that  the  assumption 
of  stationary  disturbances  gives  rise  to  two-dimen- 
sional signal  and  noise  transforms  that  cause  the 
mean  square  error  expression  to  become  independent 
of  time.  As  shown  in  f3j,  the  two-dimensional 
transforms  of  stationary  processes  do  not  exist 
unless  impulses  are  included  in  the  transforms. 

When  such  transforms  are  substituted  in  (4.10),  the. 
right  hand  side  of  this  expression  yields  an  inte- 
gral that  is  independent  of  time.  Inis  integral  is 
minimized  by  the  causal,  linear,  time-invariant 
filter  arrived  at  by  N.  Wiener.  The  sigr.31  and 
noise  terms  in  this  integral  will  be  represented  by 
the  ordinary  power  spectral  densities  instead  of 
the  average  energy  spectral  densities  which  will  be 
defined  in  Section  V. 


V.  Consideration  of  the  Integral  of  the 
Kean  Square  Estimation  Error 


Consider  equation  (4.10).  We  are  concerned 
with  selecting  an  estimator  H(ju>)  that  will  mini- 
mize this  expression  or  some  function  of  it. 
Suppose  we  address  this  problem  directly  by  trying 
to  o.inimize  (4.10)  for  each  value  of  c.  If  we 
make  the  change  of  variable 


“i  + “2 


and  hold  constant  while  integrating  with  respect 
to  p,  equation  (4.10)  becomes 


Ree(t.O  - 


to  / {i  / rss(‘Jl'P-“l) 

[i  CJoi.j)  -II  (Juj)]  [i  (J  (p-Uj ) ) -H  ( J Cp-Uj ) ) j 

+rnn  (uj  .p*v*  [H  (J“i  >]  [H  (J  {p-“i  > >] 


ddjj  J e^ptdp 


(5.1) 


We  must  now  find  H(jsu)  such  that  this  expression 
or  some  function  of 
value  of  t. 


Ree(t,t)  is  minimized  for  each 


Note  that  the  bracketed  expression  is  an  integral 


with  respect  to  and  it  will  result  in  a function 
of  p only.  If  we  denote  this  function  of  p by 
E xP).  then  (5.1)  becomes: 


J 

2it 


/ 


E2(p)eJptdp 


(5.2) 


From  (5.J)  we  Immediately  recopnizc  S<((t,t) 

£ Eie“(t)}  and  b^lp)  as  a Fourier  tiansfin  pair. 
We  denote  this  by: 


E{e2(t))  t E2(p) 

Using  this  notation,  we  also  recall: 


E2(p)  - / E(e2(t))e'Jptdt 


(5.3) 


In  particular,  for  p 
result : 


0 we  have  the  simple 


/ E{e2 (t) }dl  = E2(p) 


(5.4) 


p-0 


Note  that  (5.4)  implies  the  doubly  infinite  time 
integral  of  the  mj-an  square  estimation  error  is 
finite  provided  E‘ (p)  i is  finite.  From  (5.1) 


p-0 


e2(p)'p-o  “ H /frss(u1.-ui1)[|i(>1)-H(jcu1)|2] 

-®  ' 

•*-rnn(u1,-u.1)[|H(j.u1)|2]|du)l  (5.5) 


The  terms  rss(cui  ,-iUj ) and  rnn(uj  ,-uij)  are  called 
the  average  energy  spectral  densities  of  the 
signal  and  noise  processes.  The  right-hand  side 
of  (5.5)  bears  a striking  resemblance  to  the  mean 
square  error  expression  of  a Wiener  filter.  This 
similarity  immediately  suggests  that  one  may  be 
able  to  select  a causal,  linear,  time-invariant 
estimator  H(jui),  that  will  minimize  (5.5)  ar.d 
thereby  the  integral  of  the  mean  square  estimation 
error.  The  resulting  estimator  may  be  different 
from  the  Wiener  estimator  since  the  position  of 
the  power  spectral  density  has  been  assumed  by  the 
energy  spectral  density. 


VI.  The  Optical  Linear  Time-Invariant  Estimator 
for  a Restricted  Class  of  Nonstationary 
Processes 


From  the  previous  section  we  note  that  equa- 
tion (5.4)  and  (5.5)  yield: 


OD  OD 

/E{e2(t))dt  --L  / {rss(u.-u>) 

^|l(ju)-H(ju)  |2J+rnn(u,-u)) 

[|H(J(d)|2]  J du  (6.1) 


'in  addition,  if  we  assume  that  the  estimate 
is  unbiased,  then  we  can  use  the  terms  mean  square 
error  and  error  variance  interchangeably . 


(7.:) 


As.  Burning  tills  Integral  exists,  the  first  quest  lot; 
one  a<»ks  is  whether  or  not  there  is  an  N(Jui)  for 
which  (6.1)  Is  minimized.  This  question  Is 
answered  for  a restricted  class  of.  processes  by  the 
following  theorem,  which  applies  to  the  filtering 
case  (l(Ju)  - 1). 


lRss(tl,t2) 'dt,d‘2 


lRnn(trt2)ldtidt2  4 


Suppose  that: 

a.  The  average  energy  spectral  densities 

I*  (u.-to)  and  rnn(u),-u)  are  rational  functions 
oi  u. 

b.  The  functions  rnn(ut-iu)  have  no  repeated  poles 
in  the  u-plane. 

Then  the  choice  of  H(ju)  which  is  linear,  causal 
(al)  poles  in  the  left  half  s-plane)  and  minimizes 
the  expansion  on  the  left  side  of  (6.1)  for 
I(Ju)  * 1 (the  filtering  case)  is  given  by:. 


H(Jui)  - H(r.)| 


tryy(*)l 


f j (s)  1 


V*>  ‘ "yy<S)1+  ‘ ?ss<s>+Fnn<s> 


rss(S)  - rBs(u.-^) 


r («)  - r (u.-t) 

on  nn 


r,s(s> 

‘Vs’1' 


Fj(s)  + F2(-s) 


may  not  be  all  that  large,  but  it  is  quite  possi- 
ble that  some  of  the  processes  encountered  in 
practical  situations  may  be  approximated  b>  such 
processes. 

If  one  considers  the  nature  of  those  pro- 
cesses that  arise  from  random  initial  conditions 
applied  to  linear  dynamical  systems  for  instance, 
one  finds  autocorrelation  functions,  of  the  form 
that  satisfy  (7.1)  and  (7.2).  The  following  exam- 
ple demonstrates  the  application  of  the  theorem  to 
a three-dimensional  system  with  scalar  observa- 
tions and  white  nonstationary  disturbances. 

Example 

In  this  example.it  is  assumed  that  the  system 
is  described  by  a vector  matrix  differential  equa- 
tion. The  observation  is  a scalar  valued  function 
of  the  system's  states  with  additive  white  nonsta- 
tior.ary  noise. 

Assumptions: 

Vector  System 


.ir  1 

u iM 


Scalar  Observation 


Fj(s)  denotes  the  factored  portion  that  is 

analytic  in  the  right  one-half  s-plane 

F^C-s)  denotes  the  factored  portion  that  is 
analytic  in  the  left  cne-half  s-plane 

[ ]+  denotes  that  portion  of  [ ] having 

poles  and  zeroes  in  the  left  one-half 
s-plane  only 

[ ] denotes  that  portion  of  [ J having 

poles  and  zeroes  in  the  right  one-half 
s-plane  only. 


E{u(t)J  - E[v(t)]  - 0 

-0.l(|t.|+|t2|) 

E[u(tj)u(t2)]  • e 4(tl*t2) 


Proof 

The  proof  of  this  theorem  is  given  in 
reference  (<]  as  theorem  2. 

VII.  Application  of  the  Result 

The  result  given  by  the  theorem  indicates  that 
for  nonstationary  processes  with  two-dimensional 
Fourier  transforms  it  is  possible  to  arrive  at  an 
H(ju>)  for  which  the  lncegral  of  the  mean  square 
estimation  error  Is  minimized.  Clearly,  the  class 
of  nonstationary  processes  for  which: 


-0.5Clt.MtJ) 

E[v(tj)v(t2)]  * e d(tj-t2) 


For  thlf  system  the  transfer  function  Is  Just 


C(jui)  ’ (j-  ♦ 2)  (jw  + 3)  (ju.  ♦ 1) 


The  tve-dimenr  lnr.al  transform  of  the  signal 
autocorrelation  Is 


• . ' - N NAM  N >RES$ 

diversity  of  (California 
gineering  Systems  XLepartnent 


U.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Rose 
Bolling  AFB,  Washington,  DC  20332 


iS.  (of  this  report) 


10.  ABSTRACT  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number ) jj 

4 A technique^ is  developed  whereby  one  can  synthesize  a casual  linear  time-invariant 
estimator  that  is  optimal  for  restricted  types  of  nonstationary  processes.  The 
technique  is  applicable  to  linear,  time-invariant  systems  (driven  by  nonstation- 
ary state  noises)  for  which  scalar  observations  are  made  in  the  presence  of  add- 
itive nonstationary  noise.  Two-dimensional  Fourier  transforms  are  used  to  obtain 
an  expression  for  the  estimator's  mean  square  error.  It  is  assumed  that  it  is 
desirable  to  minimize  the  time  integral  of  this  expression.  The  calculation 
of  this  integral  results  in  an  expression  which  can  be  minimized  by  selecting 
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* 


r (u-.u..,)  - r <*., .<0  c(ju!.)  c(ju,) 

• si  .'  uu  1 i 1 * 


4(0.1)’  +(u,+u5)  j | l juj+2)  ( ju j+3)  ( juij+ 1 ) 


2 

()u2-*2)  (juj2t3)  (jw2+l) 

_ 


The  observation  noise  transform  is: 


r (u>  u ) - r (U  u.o ? 

nn  1 2 VV  1 2 / in  c\2,/  , *2 

4(0.5) 


The  average  energy  spectra  becomes: 


r (u),-ti»)  • 

(0. 1)  (uj2+(2)2)  (ui?+(3)2)  (u2+l) 


rnn<“«-»  " 03 


Using  these  sp ■ rra  the  optimal  filter,  as 
Indicated  by  the-  2 Is  just 


0.1436  (s*  + 6.071s  * 10.324) 

(s  + 3.0718)  (s2  + 3.0718s  + 2.4359) 


VIII.  Conclusions 

This  paper  has  discussed  a technique  whereby 
one  can  synthesize  a causal,  linear,  rime-invariant 
estimator  that  is  optimal  for  restricted  types  of 
nonstat ionary  processes.  The  results  are  applica- 
ble to  linear,  time-invariant  systems  (driven  by 
nonstationary  state  noise)  for  which  scalar  obser- 
vations are  made  in  the  presence  of  additive  non- 
statlonary  noise.  The  resulting  estimator  is 
similar  in  some  respects  to  the  Wiener  filter  that 
can  be  derived  under  the  assumptions  of  stationary 
state  and  observation  noise  processes.  This  simi- 
larity is  reflected  in  the  fact  that  the  power 
spectra  of  the  stationary  processes  are  replaced  by 
the  average  energy  spectra  of  the  nonstationary 
processes.  The  estimator’s  usefulness  is  limited 
by  the  requirement  that  the  observations  be  scalar, 
and  the  nonstationary  processes  have  Fourier  trans- 
formable autocorrelation  functions. 
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